Analytical estimate of effective charges at saturation in Poisson-Boltzmann cell models 
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We propose a simple approximation scheme to compute the effective charge of highly charged 
colloids (spherical or cylindrical with infinite length) . Within non- linear Poisson-Boltzmann theory, 
we start from an expression of the effective charge in the infinite dilution limit which is asymptotically 
' valid for large salt concentrations; this result is then extended to finite colloidal concentration, 

approximating the salt partitioning effect which relates the salt content in the suspension to that of 
a dializing reservoir. This leads to an analytical expression of the effective charge as a function of 
colloid volume fraction and salt concentration. These results compare favorably with the effective 
charges at saturation (i.e. in the limit of large bare charge) computed numerically following the 
standard prescription proposed by Alexander et al. within the cell model. 
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I. INTRODUCTION 



^ In colloidal suspensions, the charge ratio between the poly- ions (colloids) and the counter- ions or salt ions may be 

g 1 larger than 10*. This considerable asymmetry invalidates most of standard liquid state theories. The concept of charge 
^ ■ renormalization however allows to significantly simplify the description of such systems: the electrostatic coupling 
I— 1 1 between the poly-ions and oppositely charged micro-species induces a strong accumulation (or "condensation" ) near 
the poly-ion surface. The idea is then to consider a colloid plus its "captive" micro-ions as an entity carrying an 
^ ' effective (or renormalized) charge ZeS, a priori much smaller than the bare charge .^bare- Consequently, except in 
, the immediate vicinity of the colloids where linearization schemes fail, the interactions in the suspension are well 
' described by Debye-Hiickel-like linearized theories, provided that the bare charge is replaced by the renormalized one 
■nj" Zcs- Schematically, the linear effects of screening by an electrolyte induce a dressing of the bare Coulomb potential 
[Zbarc/f — > ■^bare exp(— Kr)/r whcrc is the Debye length, see below], while the non- linear effects of screening imply 
_ the identification Zbare exp(— Kr)/r Zcsexp{—K'r)/r (with possibly k' k [1]). 

■ Several reviews discussing the notion of charge renormalization have appeared recently [1-4]. In the colloid science 
' field this concept has been introduced by Alexander et al. [5] in the context of the Poisson-Boltzmann (PB) cell 
. model, but had been widely accepted since the fifties in the field of linear polyelectrolytes [6,7]. The definition of an 

a effective charge from the far field potential created by an isolated macro-ion in an electrolyte is unambiguous [2,8,9], 
I . at least for the simple cases of spherical or infinitely long cylindrical macro-ions we shall consider here, even if the 
'"^ ' case of less symmetric poly- ions deserves more attention (see e.g. [10]). Within a cell model, introduced to replace 
^ ' the complicated many body problem of colloids in solution by a simpler one-particle system [11-14], the definition of 
O i an effective charge is more elusive, but the celebrated proposal made by Alexander et al. amounts to find the optimal 
^ linearized PB potential matching the non-linear one at the cell boundary. The cell approach, which validity has been 
^ . assessed by several studies [16-21], appears to provide a reasonable description of solutions containing mono-valent 
counter- ions [22-24]. The effective charge is obtained by integrating the charge density deduced from the linearized 
potential over the region accessible to micro-ions, or equivalently from Gauss theorem at the colloid's surface. Note 



■ that a renormalization of Z also implies a renormalization of the screening constant k. For low poly-ion bare charges 
■^baro, linearizing PB equation is a valid approximation so that Z^s ~ 2^baro, whereas in the opposite limit of high 
bare charges, Z^s was found numerically to saturate to a value independent of ^baic: ^cft — ^sat- The saturation 
value Zsnt depends on the geometry of the colloid, the temperature or the quantity of added-salt. Unfortunately, no 
analytical prediction is available for these dependences. 

In this letter, we propose an approximate analytical expression for the effective charge of spherical or rod-like 
macro-ions, as a function of macro-ion density or salt content. Our approach starts with an estimate of the saturation 
value of the effective charge in the infinite dilution limit, as deduced from recent analytical results [25]. In the case 
of finite colloidal dilution, the osmotic equilibrium of the suspension with a salt reservoir is modeled using a Donnan 
equilibrium approximation. This allows to derive a simple polynomial equation of degree 4 fulfilled by the saturation 
value Zsat of the effective charge Zcs, in the case of a symmetric 1:1 electrolyte. In the limiting case of no added salt. 
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the above equation is easily solved analytically and provides results in quantitative agreement with the saturation 
effective charges following Alexander's prescription. 



II. GENERAL FRAMEWORK AND METHOD 



We consider first the situation of an isolated macro-ion of given surface charge density in a electrolyte of bulk 
density no (no confinement). The solvent is considered as a medium of uniform dielectric (CGS) permittivity e. 
Within Poisson-Boltzmann theory, the electrostatic potential, when assumed to vanish far from macro-ion obeys the 
equation 

V^(/) = Kosinh^, (1) 

where the screening factor kq is defined as Kq = SttIbtio and the Bjerrum length quantifies the strength of electrostatic 
coupling: £b = e^/{skT) (e > denotes the elementary charge and kT is the thermal energy). A complete asymptotic 
solution of Eq. (1) has been obtained recently by Shkel et al. [25] in spherical and cylindrical geometry. Prom the 
far-field behaviour of the corresponding solutions, one obtains after some algebra the effective charges at saturation 

Zsat — = 4:Koa + 6 + 1 — — I for spheres (2) 

a \KoaJ 

Asat^B = 2Koa + i + ( — ) for cyhnders. (3) 
2 \KoaJ 

In these equations, a denotes the radius of the macro-ion under consideration, Ze the total charge in the case of 
spheres and Ae the line charge density in the case of cylinders (with infinite length). Expressions (2) and (3) are the 
exact expansions of the saturation charges in the limit of large kqu but are in practice accurate as soon as kqu > 1 
(not shown). 

Our goal is to translate these relations into expressions that would approximate the effective charges in confined 
geometry, where the macro- ion is enclosed in a cell [5]. To this aim, we first find an approximation for the relevant 
screening factor before inserting it into (2) and (3) making the substitution Ki, kq- The implicit assumption is 
that the mean salt density in the cell is related to the effective charge in a similar manner as in the infinite dilution 
limit. 

In confined geometry, PB equation still takes the form (1), where kq = {SnlBnoY^'^ should now be considered as the 
inverse screening length in a (neutral) salt reservoir in osmotic equilibrium with the solution, through a membrane 
permeable to micro-species (no is thus now the salt density in the reservoir). It may be shown that the relevant 
screening factor in the cell is related to the micro-ions density at the cell boundary [8] : 

kI = 4TTiB [p+ (RWS) + P~ {RWS)] ■ (4) 

In this equation Rws is the radius of the Wigner-Seitz (WS) confining cell. We have recently proposed an efficient 

prescription to compute without solving the complicated non- linear problem [8] . It is however impossible to deduce 
an analytical expression of from this approach and we resort to the following approximation. We assume that the 
micro-ion densities are slowly varying in the WS cell so that the mean densities provide a reasonable estimation 
of the boundary densities p^{Rws)- We thus write 

kI = 47r£B(n+ +n-) = A-KiB{2n+ + Z^^tp), (5) 

where the last equality follows from the electroneutrality constraint (p is the density of colloids, assumed positively 
charged without loss of generality). Note that at the level of a linear theory, the effective charge Zgat and not the 
bare one Z^aie enters this expression. 

We now need to relate the mean densities to no, the concentration in the reservoir (the so called Donnan effect, 
see e.g. [26,14]). Chemical equilibrium imposes p'^p~ = Uq at any point in the cell. We again assume this relation to 
hold for the mean densities, so that n+n" = Uq. This leads to 

Kt = 4 + (47r£B^satp)^ (6) 

Up to now the reasoning is quite general and independent of the geometry. For spheres with radius a and packing 
fraction t] = Anpa^/S = (a/Rws)^, and for rods with packing fraction r] = (a/Rws)^, we obtain 

(«;*a)^ = (koo)* + (3r?.^sat^B/a)^ (spheres) and («;*a)* = (koo)* -|- (4ryAsat^B)^ (cylinders) (7) 
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Supplementing these equations with the k dependence of Zgat obtained in Eq. (2) for spheres and (3) for cyUnders, 
leads to an equation of degree 4 satisfied by n^, (or equivalently ^sat)- 

This equation can be solved analytically. We give the general solution in appendix A. However, without added salt 
(formally kq = 0), the solutions of these equations take a particularly simple form: 

^sat ^ = 6 + 24r7 + 12 + V (8) 
Asat^B = |+8r? + 2V6r? + 16r;2 (9) 

We compare on Figs. (1) and (2), the results obtained following this route with the full PB estimate following 
Alexander's prescription. The latter involves a numerical resolution of the PB equation. Both the no-salt and the 
finite ionic strength case are considered. The agreement is seen to be quite good in view of the minimum number of 
ingredients involved in the present approach. 



III. CONCLUSION 



In this contribution, we propose a simple approximation scheme to compute the saturation value of the effective 

charge in concentrated suspensions of highly charged colloidal particles. Our estimate starts with an asymptotic 
expression of the effective charge in the infinite dilution limit, obtained from recent analytical results [25]. The case 
of finite colloidal dilution is described using a Donnan equilibrium approximation for the osmotic equilibrium of the 
suspension with a salt reservoir. This calculation leads to an analytical estimate of Zg^^t, the saturation value of the 
effective charge, as a function of the density of colloids and salt concentration. Our starting points [Eqs. (2) and (3)] 
neglect contributions of order (koo)"^ and become inaccurate for na < 1. This implies that our effective charges at 
finite density becomes less reliable in the salt free case for small packing fractions rj. However, this is seen to occur 
for very small rj only (see e.g. Fig. 1). 

These results are obtained at the level of the mean-field Poisson-Boltzmann theory. We note however that the 
existence of a saturation value of the effective charge, .^sat) independent of the bare charge, is indeed confirmed 
in more refined approaches in the colloidal limit a ^ Ib, with a the colloid size and £b the Bjerrum length (see 
e.g. Groot [22] using the primitive model). These results show moreover that PB theory becomes successful in the 
afore-mentioned colloidal limit, a 3> [22,20]. Eventually the saturation picture within PB theory is in quantitative 
agreement with experimental data for the osmotic pressure [27,28,8]. 

It is eventually instructive to reconsider the results reported by Alexander et al. in their original paper [5]. At a 
packing fraction rj = 0.125, they find numerically a saturation value for ZeS of the order of 15, in a/Is units (spherical 
coUoid). On the other hand our expression (8) for the same r] gives Z^^ilsja = 9 -|- 2>\fh ~ 15.7, i.e. very close to the 
value found in [5]. Note that this number of 15 has subsequently been often quoted in the literature as a "standard" 
value of the effective charge. The saturation value of the effective charge however crucially depends on the colloid 
volume fraction, and the full density dependence of Zgat should be taken into account in finite concentration cases. 
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FIG. 1. Effective charge at saturation Zsnt (in units of radius over Bjcrrum length £b), as a function of packing fraction 
for a spherical macro-ion of radius a enclosed in a concentric spherical cell of radius arf^^^'^ . The situation is without added 
salt. The analytical expression (8) is shown by the continuous curve while the non-linear Poisson-Boltzmann values computed 
numerically following Alexander et al. [5] are represented with circles. Inset: same with added salt for Koa = 2.6. 
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FIG. 2. Effective line charge density at saturation Asat (in units of as a function of packing fraction for a cylindrical 

macro- ion of radius a enclosed in a concentric cylindrical cell of radius arj~^^^ (no added salt). The analytical expression (9) 
is shown by the continuous curve while the non-linear Poisson-Boltzmann values computed numerically following Alexander et 
al. [5] are represented with circles. The inset shows the same quantity with added electrolyte (kqo = 3.0). 



APPENDIX A: ANALYTICAL EXPRESSION FOR THE EFFECTIVE CHARGE 



Wo give hero the explicit expression for the effective charge as a function of the colloid packing fraction r] and salt 
concentration in the reservoir, no = Kq/Stt^b. 

Inserting Eqs. (2) or (3) into Eq. (7), one gets the quartic equation for X = K-^a : 

X'' = X^ + Tf{a + (3Xf (Al) 

with Xq = Koa. Values for a,/3 are respectively a = 12 and /3 = 18 for spheres, and a = 8 and /? = 6 for cylinders. 
Let us introduce the following quantities : 



A = -12X0^ - 12(3^ rf + a*??^ 

r = 108a2/3\* - 2aS' - T^ocWiXt + PW) 

$ = -72X^a^rf + 36a^p'^r]'^ - 2a^rf 



The solution X = Ki,a is then found as 

1 /2a2^2 ~ I /4q,2 2 2 

K^a = --1/ + @+ - { Oh , 

2V 3 2V 3 



e 



1/2 

(A3) 



The effective charge is then obtained by replacing this value for in Eq. (2) for spherical macro ions, and in Eq. 
(3) for rod-like macro ions. 
In the no-added salt case, these expressions reduce to Eqs. (8) and (9), as indicated in the text. 



5 



